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Phase Separation and Coarsening in Electrostatically Driven Granular Media 
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A continuum model for the phase separation and coarsening, observed in electrostatically driven 
granular media, is formulated in terms of a Ginzburg-Landau equation subject to conservation of 
the total number of grains. In the regime of well- developed clusters, the continuum model is used 
to derive "sharp- interface" equations that govern the dynamics of the interphase boundary. The 
model captures the essential physics of this system. 
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Despite extensive work in the last two decades, the 
physics of granular flow is still poorly understood, es- 
pecially in the limit of strongly inelastic granular colli- 
sions when no first-principles hydrodynamic description 
is available jj],[§. Additional complications include con- 
tact interactions which become dominant when the grain 
size goes below 0.1 mm. As particles may acquire an 
electric charge, a new type of dynamics appears which 
is governed by the interplay between long-range electro- 
magnetic and short-range contact forces. In this regime, 
electrostatic excitation of granular media becomes possi- 
ble. It offers an opportunity for an investigation of gran- 
ular systems consisting of small particles. 

Recently, off-equilibrium phase separation and coars- 
ening in electrostatically driven granular submonolayers 
was reported ||,[|. The inset of Fig. 1 shows a schematic 
of the setting. Conducting particles are placed between 
the plates of a plane capacitor which is energized by a 
constant (DC) or alternating (AC) electric field E. When 
particles are in contact with the capacitor plate, they ac- 
quire an electric charge. If the electric field exceeds a 
critical value, the resulting (upward) electric force over- 
comes the gravity force mg (g is gravity acceleration and 
m is mass of particle) and pushes the particles upward. 
When a particle hits the upper plate, it gets the opposite 
charge and falls back. Changing the frequency uo of the 
AC field E(t) = Eq cos(o;t), one can control the particle 
elevation and avoid the collisions with the upper plate. 

It was found in Ref. || that the particles remain im- 
mobile on the bottom plate at E < E\ (the precipitate 
state). If the electric field is larger than a second thresh- 
old value, E 2 > Ei, the system is in a gas-like state. 
This second field E 2 is 50%-70% larger than Ei. It was 
found that, upon decreasing the field below E 2 (in the 
interval E\ < E < E 2 ) there is nucleation of precipi- 
tate, and small densely packed clusters form and start 
to grow. The clusters then exhibit Ostwald ripening ||: 
smaller clusters shrink, while larger clusters grow at their 
expense. Molecular dynamics simulations showed quali- 
tative agreement with experimental results on phase sep- 
aration and coarsening 0. 



In this Letter we develop a continuum description 
of these phenomena. The model captures the essential 
phenomenology, correctly reproducing different possible 
morphologies of coarsening and providing a correct quan- 
titative description of its dynamic scaling behavior. 

Let the granulate consist of N ^> 1 identical spherical 
particles with mass m and radius cr (which is small com- 
pared with the plate spacing h). The continuum version 
of the model employs two fields: the densities of the pre- 
cipitate, n(r, t) and gas phase n g (r, t), where r = (x,y) 
and x and y are the horizontal coordinates. For con- 
creteness, we will consider a low- frequency electric field 
(so that the gas density is almost independent of the ver- 
tical coordinate) and measure n g in cm -2 . 

The phase separation is caused by electrostatic screen- 
ing ||: a decrease in the vertical electric force F, exerted 
on a grain in contact with the bottom plate, caused by 
the presence of other grains. Thus, the F(n) dependence 
is an important element of our model. This dependence 
can be easily found in the dilute limit na 2 <C 1. Con- 
sider two grains lying on the bottom plate the distance 
L ^> a apart. The second grain and its mirror image in 
the plane form a dipole with the dipole moment of or- 
der of a 3 E. This dipole produces an electric field at the 
location of the first grain: 5E ~ —(<j/L) 3 E. Therefore, 
the net force acting on the first grain is reduced and be- 
comes Fq [l — tv(a/L) 3 )] . If there are many grains lying 
far apart, we can simply sum up these vertical forces and 
obtain the net force 



F(n)=F (l-k 1 a 3 n 3/2 ). 



(1) 



Here Fq = 1.36 . . .a 2 E 2 is the vertical force exerted on 
a single grain in contact with the bottom plate ||, and 
ft, k\ = 0(1) are numerical factors. 

Although no analytic expression for F{n) is presently 
available for intermediate and large n, it is clear that 
F(n) should decrease with n [||. A decreasing F(n) de- 
pendence leads, at intermediate values of E : to a segrega- 
tion instability and bistability. Indeed, the density value 
n = n* such that F(n*) = mg is in an unstable equi- 
librium. For n < n*, F(n) exceeds mg so the particles 



"evaporate" until the "empty state" n = is reached. If 
n > n*, F(n) < mg and the particles remain immobile. 
However, as gas particles hit these regions and (because 
of strong inelasticity of collisions) get attached to them, 
n grows until the densely packed state n = n c ~ o~~ 2 
is reached. This simple argument ignores conservation of 
the total number of particles and cluster edge effects that 
we account for in the following. 

The vertical force exerted on a particle decreases with 
an increase of the number of neighboring particles [0: 
the force exerted on a particle located at the cluster edge, 
F e , is larger than the force on a particle in the bulk of 
the cluster F5, but smaller than the force on a single 
particle Fq. In addition, in a coarse-grained description 
(which should be valid for clusters with many particles) 
F e should depend on the local curvature of the cluster 
edge. Although no general relation for F e is presently 
available, a clear signature of these cluster edge effects 
appears already in the dilute limit where quantitative 
relations can be obtained. First, for an inhomogeneous 
density distribution, the density n(r) in Eq. (|l|) should be 
replaced by a locally averaged density n(r) , the averaging 
being performed over a region which size is of the order of 
the interparticle distance. For a weakly inhomogeneous 
coverage, one can expand n(r) around a point 1*0 up to 
(r — ro) 2 . Then, averaging the result over a circle of 
radius ~ n(r ) -1 / 2 , we obtain n(r) = n(r) + 0(V 2 n/n). 
Substituting it in Eq. (|l|), we arrive at 

F(n, V 2 n) = F (l - jfei a 3 n 3 / 2 - k 2 a 3 n~^ 2 \/ 2 n) , (2) 

where k 2 = 0(1) is another numerical factor. The cluster 
edge effects are described by the V 2 -term. 

Now we start formulating our model. The precipitate 
number density n(r, t) will serve as the order parame- 
ter, the two phases corresponding to n = and n = n c 
||. In its turn, the gas density n g (t) plays the role of 
the "mean field". (We assume that the density relax- 
ation in the gas phase is fast compared to the cluster-gas 
exchange dynamics. Therefore, n g is approximately con- 
stant in space and depends only on time.) First, consider 
a uniform precipitate in the dilute limit: n + n g <C n c . 
The precipitate will "evaporate" if Fq > mg. On the 
contrary, the gas will precipitate if Fo < mg. The typical 
time scales of each of these two processes are of order of 
r = min [(g/h) 1 / 2 , a; -1 ]. A simple dynamic equation 



^ = ^[Cn g 0(T)+n0(-T)}, 



(3) 



subject to conservation law n + n g = N/L 2 = const, ac- 
counts for these facts. Here T = 1 — Fo/mg and 0(. . .) is 
the #- function. The numerical factor C = 0(1) accounts 
for a possible difference between the time scales of "evap- 
oration" and precipitation. The factor T in front of the 
brackets of Eq. (^) takes care, in the simplest possible 
way, of the continuity of transition between the regimes 
of Fq > mg and Fq < mg. 



For higher precipitate densities one needs to account 
for additional effects. First, because of the screening, the 
lifting force decreases with n/n c . For simplicity, we as- 
sume a linear dependence Fq(1 — bn/n c ), for all densities 
< n/n c < 1. Here b is a numerical constant (using the 
estimate of F(n = n c ) from Ref. Q, we have b c^ 8/9). 
Second, we account for the inhomogeneity of the precip- 
itate by adopting the V 2 n term from Eq. (||), with a 
diffusion coefficient D ~ a 2 jr. Third, we introduce an 
additional factor (n c — n)/n c in the first term in the r.h.s. 
of Eq. (ph. It accounts, in a simple way, for a slowdown 
of precipitation from the gas phase: at finite n/n c a part 
of the bottom plate is already occupied by grains. 

In this way one arrives at a scalar Ginzburg-Landau 
equation (GLE). We will be using this equation in the 
phase coexistence regime mg < Fq < mg/(l — b). (A 
more detailed condition will be derived below.) The 
equation can be written in a scaled form: 



dn 
~dt'~ 

where 

(j) = (n — n*) x 



cj){n,n g ,n^) + V 2 n , 



/ n, 

I Cn g (l 



if < n < n* 
n), if n* < n < 1 . 



(4) 



(5) 



Here n* = (1/6) (1 — mg/Fo). The coordinates (and 
the system size L, see below) are scaled by S = 
( J Dr/A) 1 / 2 /n c , the time is scaled by rmg/bFo. The pre- 
cipitate and gas densities (and n*) are scaled by n c . Fi- 
nally, A = bF^/mgn 2 . 

The dynamics are constrained by the conservation of 
the total number of particles. In the scaled units 
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n(x,y, t) dx dy + n g (t) = e , 



(6) 



where e is the (constant) area fraction of the granulate. 
At fixed n g and < n* < 1, the function <j)(n,n g ,n*) 
describes bistability: it has two stable zeros: at n = and 
n = 1, and an unstable zero at n = n*. Eqs. (^-(Ig) serve 
as the continuum version of our model. Similar equations 
(with different forms of function (j) and conservation laws) 
have appeared in other contexts [0-0 . The dynamics of 
a globally conserved bistable system involve three char- 
acteristic time scales. Phase separation occurs on the 
fastest time scale to, determined by the properties of the 
function 0; it is independent of the characteristic cluster 
size /. On a longer time scale t\ ex /, the time-dependent 
gas density n g (t) becomes close to special "equilibrium" 
value n eq such that the area rule 



(n,n g ,n*) dn = 



approximately holds. For the present model, a straight- 
forward calculation yields 



-.eg 



C(l-n*) 3 



(7) 



When n g = n e g q , a planar interface of the precipitate is in 
a dynamic equilibrium: it neither advances nor retreats 
if fluctuations are neglected. For curved interfaces there 
is an additional, slower coarsening stage of the dynamics, 
governed by the cluster edge curvature. Its characteristic 
time t 2 ocl 2 §. 

A detailed analysis of this system is possible in the 
asymptotic sharp-interface limit. The sharp-interface 
equations become valid towards the end of the ti-stage, 
when the precipitate in the clusters is already densely 
packed (n = 1), and the gas density n g is close to n e g q . In 
this limit, the conservation law (pi) reads 



L 2 



(8) 



where A(t) is the total area of the precipitate. Demand- 
ing A > and using Eq. (0), we obtain a condition for 
the two-phase coexistence: 



e > 



C(l-n*) 3 * 
Using our model expression for n*, we can write 



mg 



Ml 



(9) 



(10) 



One should compare the more restrictive condition (|io| ) 
with the "naive" condition 1 < Fo/mg < (1 — ft) -1 . For 
example, for e = 0.25, C = 1, and b = 8/9 the "naive" 
condition gives 1 < Fo/mg < 9, whereas Eq. (|l(j) gives 
1 < Fo/mg < 1.5. The latter condition agrees better 
with experiment Q. Even better agreement is achieved 
if one chooses C — 10. 

The normal component of the interface speed is |1Q] : 



vC 



(1-n,) 3 

3/2 



K 



? )-/c, 



(11) 



where /C is the local curvature of the interface and v = 
5(9 + 2>/3)/138 « 0.45. 

Given the initial location of all interfaces, Eqs. (g) 
and (O) provide a proper description of the late-time 
coarsening dynamics. Furthermore, these equations can 
be mapped into equations of interface- controlled trans- 
port which appeared in other contexts. This enables one 
to readily present a number of important results for sev- 
eral coarsening configurations: (1) A planar or circular 
interfaces are stable with respect to small modulations 
P,PT[; (2) The radius of a single cluster which can be in 
a stable dynamic equilibrium with the gas phase has a 
non-zero lower bound. This lower bound scales with the 
system size (like L 2 / 3 ), but is independent of N |]9|,pJ; 



(3) Multiple clusters exhibit Ostwald ripening: j7|-p|,pT[; 

(4) If there are many clusters, their size distribution ex- 
hibits dynamic scaling. The number of clusters decays 
with time like t~ 1 , while the average cluster radius grows 
like t 1 / 2 P ]ll| , |i2| ] ; (5) Depending on the initial parame- 
ters, a "hole" (empty region) inside a cluster will either 
shrink to zero, or expand and come out of the cluster [ [l3| . 
It should be stressed that these results are insensitive to 
the exact form of the bistable function cf). 

Properties 1, 3 and 4 were already observed in experi- 
ment P] . We checked predictions 1 , 2 and 5 by perform- 
ing additional experiments with cells similar to those de- 
scribed in Ref. (3). The experimental results fully agree 
with these predictions. 
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FIG. 1. The gas phase density n g (normalized to n c ) vs. 
time. Inset: schematic of the experimental setting H. 

Though suitable for theoretical analysis, sharp- 
interface equations are less convenient for numerical sim- 
ulations. We solved Eqs. (§)-@ numerically in peri- 
odic boundary conditions. Figure 2 shows a typical sim- 
ulated evolution of the precipitated phase. The initial 
conditions describe a situation when most of the par- 
ticles are in the low-density precipitate (that is, in the 
unstable region). This corresponds to an up-quench of 
the electric field from E < E\ to the coexistence region 
E\ < E < E2. At small times, "holes" in the precipitate 
develop and grow. At later times multiple clusters form 
and exhibit Ostwald ripening, in agreement with experi- 
ment. Figure 1 shows the time history of the gas density 
rig (normalized to n c ). At late times, n g approaches the 
value corresponding to the coexistence of a single cluster 
and gas phase |9|,|ll||. For the set of parameters chosen 
for this simulation, the area rule value n e g q = 0.05. The 
fine structure of the n g (t) dependence (small peaks) cor- 
responds to the disappearance of clusters Qjlj . In ex- 
periment, the current through the cell is carried by the 
grains belonging to the gas phase, so it should be pro- 
portional to n g . Therefore, the "mean field" n g (t) is a 
directly observable quantity. 

The "up-quench" type of morphology at small times 



shown in Fig. 2 was not reported in Ref. ||, so we per- 
formed special experiments with electrostatic cells sim- 
ilar to those used in Ref. |J. Starting from the gas 
phase (E > E2) and making a down-quench to E < E\, 
we first prepared a nearly uniform layer of precipitate. 
Then, after an up-quench into the coexistence region 
Ei < E < E2 , we observed the same type of morphology 
as shown in Fig. 2. 
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FIG. 2. Simulated dynamics of the precipitate for t = 100 
(a), 160 (b), 300 (c) 600 (d) 4000 (e) and 6000 (f) scaled time 
units. White color corresponds to densely packed clusters 
(n ~ ra c ), black to empty regions. The system dimensions (in 
the units of 5) are 500 x 500. The electric field E corresponds 
to n* = 0.2. 

It should be noted that when F is only slightly larger 
than mg, the interface width (which, at n* <C 1, is 
proportional to n* ' ) becomes comparable to the clus- 
ter size and/or intercluster distance. In this regime the 
sharp-interface theory is invalid. Still, the continuum 
theory should work. 

The present model assumes a low-frequency electric 
field. In this case the gas particles typically perform an 
appreciable bounce motion, and a separate, mean-field 
treatment of the gas phase is legitimate. For very high 
frequencies, the particle motion becomes effectively re- 
stricted to the bottom plate. The corresponding coars- 
ening dynamics can be quite different. One can expect 
that, at the high frequencies, the 1/2 growth exponent 
(observed in experiment || and predicted by our model) 
will cross over to the 1/3 growth exponent, typical for 
/oca%-conserved systems Q. 

It is interesting to compare the phase separation prop- 
erties of electrostatically driven granulates with those 
vibrated in the vertical direction mechanically [[15|. 
Though these two types of systems strongly differ in de- 
tails of particle interactions and motion, they have very 
similar phase diagrams and are strikingly similar in their 
phase transition morphologies. This gives an indication 
that a bistable Ginzburg-Landau equation subject to con- 
servation of the number of particles can be relevant for 



the mechanically vibrated systems as well. 

In summary, we have formulated a phenomenologi- 
cal theory of the dynamics of off-equilibrium phase sep- 
aration in electrostatically driven conducting monodis- 
perse particles. The theory includes both continuum, 
and sharp-interface formulations. It captures the essen- 
tial physics of this system. Several new predictions of 
our model have been verified in experiment. A better 
quantitative understanding of this system requires ad- 
ditional quantitative experiments in different regimes of 
phase ordering. A comparison of the model with the new 
experiments will enable one to determine the unknown 
numerical coefficients of the model. These coefficients 
control the precise phase diagram, the interface velocity, 
the amplitudes of the dynamic scaling laws in the Ost- 
wald ripening regime, etc. However, already at this stage 
we see a strong evidence in favor of quantitative relevance 
of the globally constrained GLE to this system. 
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